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Abstract 

In this paper, we study gradient decay estimates for solutions to the multi¬ 
dimensional Fisher-KPP equation with fractional diffusion. It is known that 
this equation exhibits exponentially advancing level sets with strong qual¬ 
itative upper and lower bounds on the solution. However, little has been 
shown concerning the gradient of the solution. We prove that, under mild 
conditions on the initial data, the first and second derivatives of the solution 
obey a comparative exponential decay in time. We then use this estimate to 
prove a symmetrization result, which shows that the reaction front flattens 
and quantifiably circularizes, losing its initial structure. 

Keywords: Fisher-KPP equations, fractional diffusion, gradient decay esti¬ 
mates, traveling fronts. 


1 Introduction 

The goal of this paper is to understand a strong symmetrisation phenomenon, ob¬ 
served in [7], for the level sets of react ion-diffusion equations of Fisher-KPP type 
with fractional diffusion. The model under consideration is 

Am + A“m = /(m) (1.1) 

on M'^. Here, a G (0, 2) and A“ is the fractional Laplacian of order a/2\ 

(-A)“A„(a;) = = lirncjv,a [ 

Ay\>e \y\ 

The nonlinear term / is assumed to be of KPP-type [in]: /(O) = /(I) = 0, f{y) > 0 
for y G (0,1), and f'{y) < f{y)/y- Not at ion-wise, let k = /'(O) and X = K/{d + a). 
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Under these assumptions, and for a compactly supported initial datum u{., 0) := Uq, 
it is well-known that the level sets of u will spread exponentially fast in time (Cabre- 
Roquejoffre [5]): 

Theorem 1.1 Under the above assumptions, 

• For all c> X, we have lim sup u{x,t) = 0. 


• For all c < X, we have lim inf u(x,t) = 1. 

t—^+00 

Thus, when renormalized by the exponential, the level sets of u asymptotically look 
round. It is then natural to ask whether this property holds in a more precise 
fashion, and a first answer is that given in Cabre-Coulon-Roquejoffre [6]: for a given 
h G (0,1) we have 

{u = h}F < 1x1 < 

for some constant C > 0. The next question is whether one can make this constant 
more precise. Given the rapid growth of the level sets, one could expect the possibil¬ 
ity for very erratic behavior. To confirm this, the case was investigated numerically 
by A.-C. Coulon in her PhD thesis. We reproduce here a sample of her simulations. 
The initial datum (Fig. 1) is pictured below. 


Figure 1: Level sets of uq 

and the time evolution is shown in Fig. 2. Surprisingly, this strongly suggests 





Figure 2: Level sets of m: a = 1 (left) and a = 1.6 (right) 
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asymptotic symmetrization. So we may ask whether a very general theorem of 
Jones [1] applies. It concerns the solutions of (1.1) with a = 2: 


dtu - Au = f{u) (1.2) 

with any nonlinearity / - not limited to KPP. The theorem asserts that, if mq is 
compactly supported, then, for any later time, and for any regular value h of m, and 
for any x on the level set {u{.,t) = h}, the normal line Aj.* to the level set passing 
through X intersects the convex hull of suppwo. If the level sets of u expand, this 
has very important consequences: the level sets of u symmetrize asymptotically, and 
have bounded oscillation. This is a very remarkable result, one reason being that it 
does not depend on the precise expansion rate of the level sets of u. In particular, 
if the nonlinearity / is of KPP type, then 

{u = h} = {|a;| = 2^/Kt - ^ ^ ^ Logt + 0(1)}. (1.3) 

The term 2^/Kt is due to Aronson-Weinberger [1], and the logarithmic correction to 
Bramson [ 1 ] (d = 1), and Gartner [8] (d > 1). 

Let us briefly describe a proof, due to H. Berestycki [ 2 ] , of Jones’s theorem: assume 
the contrary, then there is a hyperplane H containing A^. ^ and such that the convex 
hull of supptio lies strictly in H~, the lower half space bounded by H. If u is the 
reflection of u{.,t) about id, and v = u — ii, then v satishes a linear equation and 
is positive in id“, since it is positive at f = 0 and vanishes on H. At time t, any 
derivative of n in a direction normal to H is nonzero, contradicting that A^. ^ C id. 
We immediately see that this argument cannot be applied to our case, simply be¬ 
cause V would need to be nonnegative in instead of id“, which is impossible. 


We are going to show that, nevertheless, a rather strong form of symmetrization 
occurs. The ingredient is the following gradient estimate for u; we believe that it is 
of general interest. 


Theorem 1.2 Assume u ( a ;, 0) := uq{x) to he continuous, nonnegative and nonzero, 
and Uo{x) = as |a;| —)■ J-oo, for some e > 0. Then we have, for a solution 

u{x,t) of (1.1), a universal constant C, and a 6 > 0 (depending on X and a): 

\'Vu{x,t)\ < Ce~^^u{x,t), (1.4) 

and 

\'V‘^u{x,t)\ < Ce~^^u{x,t). (1.5) 


Theorem allows us to prove an analogous estimate for the fractional Laplacian. 
This permits us in turn to reduce the problem to the ODE ft = (1 -|- 0{e~^^))f{u), 
which is much simpler to study. The main result of the paper is 


Theorem 1.3 Letuo he as in Theorem 1.2. For every h G (0,1), there is a constant 
gh{uo) > 0 and J > 0 such that 


{u = h} = {|x| = qhiuo)e^\l 0{e '^*))} 


( 1 . 6 ) 
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This provides an explanation to the observed behavior of u. However we note that 
onr results could be improved in two ways. 

• The assumptions on uq seem slightly non optimal: indeed, it would certainly 

suffice to have uo{x) = with e > 0. Our theorem would remain 

valid under that assumption, at the cost of heavier computations. However, 
when uq decays like (or slower than) different phenomena occur, as was 

observed in [7]. 

• We do not go as far as proving a Jones type theorem. Indeed what we obtain is 
a precise estimate of the normalized level sets instead of the true level sets. In 
other words, lower order (but still exponentially growing) terms may prevent 
a full symmetrization. 

These issues will be investigated in a future paper. 


In Section 2, we gather some known (but useful) facts that will be used throughout 
the proofs of Theorems M and |1.3[ Section 3 is devoted to the proof of Theorem 


1.2[ while Section 4 is devoted to estimating the fractional Laplacian. Theorem 1.3 


is then proved in Section 5. 


2 Preliminary material 


The gradient estimate ( |1.4[ ) will be obtained by examining its representation in three 
different ranges, which is reflected by the collection of results below, that we recall 


for the reader’s convenience. The proof for the estimate on second derivatives (1.5) 


follows a similar approach, but will itself make use of (1.4). 


2.1 Invariant coordinates 

The starting point is the following estimate, proved in [B] 
Theorem 2.1 We have, for a universal C > 0: 

c-i .... C 


1 + e 


—/ct 


X 


d-\-a. 


< u{x, t) < 


1 + e 


— K.t 


X 


d-\-a ' 


This motivates the introduction of the invariant coordinates f = xe 

dtu - Xf ■ Vgu + - f{u) = 0. 

For any hxed coordinate Xi, set 0 = d^iU. We then have 

^40 += f{u)4>. 

Letting v = d^M = yields 

dtv — Xf ■ Vv + e“"'^*A"n — f{u)v — Xv = 0. 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


(2.4) 
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For another fixed coordinate Xj, set ip = dx 4> = and p = u, with v = 

This then gives us 

dtpJ + + f'{u)pj, (2.5) 


along with the associated equation in exponential coordinates: for V = 
have 


dtV -X^-VV + - f'{u)V - f"{u)vv - 2A1/ = 0. 


d^.d^^u, we 


( 2 . 6 ) 


A consequence of Theoreni |1.3| is the large time convergence of the initially compactly 
supported (or initially exponentially decreasing) solutions of (2.2) to a steady profile. 
Notice indeed that the steady equation 


- Ae ■ = f{u) (2.7) 

has a unique one-parameter family of radial solutions (mt(|'C|))t>o- Then, assuming 
Theorem 11.31 

Theorem 2.2 There exists a Too{uo) > 0 such that 

lim u{i,t) = n^^(„o)(0, 

C—)- + 00 

uniformly on compact subsets o/M'^. 


2.2 Heat kernel 


Let pa{x,t) be the heat kernel of A", in other words the solution of pt + A"p = 0 
having, as initial datum, the Dirac mass at 0. 


Proposition 2.3 There is C > Q such that, for large t and x, we have Pa{x,t) = 
t~'^l°‘Pa{xjt'^^°‘), with 

c-^ C 


ICI 


d+a+1 


< |Vp„(C)l < 


ICI 


d+a+l’ 


( 2 . 8 ) 


and there zs c„ > 0, 5 > 0 such that 


<PaiO 


ICI 


d-\-a-\-5 


ICI 


d+Q — 


c 


ICI 


d+a+S ' 


(2.9) 


It is important to note that the estimates (2.8) and (2.9) are intended for ( bounded 
away from 0; pa and its derivatives are bounded functions (in x) for any f > 0. A 
standard way to prove the above proposition is to write 


pa{,x,t) = ^ (^e , 

and to evaluate the inverse Fourier transform with the aid of Polya integrals; see for 
instance m- 
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2.3 Comparison principles 

We will also require the following easy extension of the Maximum Principle: 


Theorem 2.4 (Selective Comparison Principle) Let Q{t) be a time-dependent 
family of compact domains in (in the variable) with smooth boundaries and 
continuous time dependence; that is, {(^, t) | ^ G is an open set in x IR+. 

Let v{f,,t) be the solution to (2.4) and let w{f,,t) be a positive function such that 
w > |n| at t = 0. If, for all t > 0 we have w > |n| on the closure of the complement 
of Q{t) and 

dtw - AC Vw + - f{u)w - Xw > 0 


everywhere inside LL{t), then w{f,,t) > |n(^,t)| on all o/M'^ x 


Proof. We first show that w > v hj looking at the equation for q = w — v on Q{t) 
which, by assumption, starts positive at t = 0. Observe that, if q achieves a global 
minimum value of 0 at time f and location ^ G Ll{f), then we easily have dtq{f,, t) > 0 
(all other terms on the left are nonpositive); this, however, crucially requires that 
w > V outside So, by continuity in time, q can never become negative inside 

fl{t). To complete the proof, we need to show that w > —v. But this is done in 
precisely the same manner as before, now examining the equation for q = w -\-v and 
showing that it too remains positive for all time. • 


Also note that the above result (with a minor modihcation) also holds for a solution 
V of (2.6). That is, if W is a positive function such that W > \V\ on the closure of 
the complement of Q{t) for all t > 0 and also W > \V\ on ^2(0), and 

dtW - Ae ■ VW + - f{u)W - 2AW > \ f\u)vv\ 


everywhere inside Lift), then W{f,,t) > |1A(^,t)| on all of x IR+. The proof is 
exactly the same since the equations for W — V and W -\-V have nonnegative inho¬ 
mogeneities on the right side of the inequality. 


Finally let us mention the following version of Kato’s inequality and its well-known 
consequence for the Fisher-KPP equation. 

Proposition 2.5 If u{x) is smooth, then 

A“|M(a:)| < sgn(M)A"n(a;), 

in the distributional sense (and in the classical sense if a <1). 

Proof. Recall the elementary inequality |a| — |6| > sgn(6)(a — h). This implies, for 
all {x, h) G X 

\u{x h)\ — |n(a;)| > sgn(M(a;))(n(a; -\- h) — u{x)). 

Integrating both sides in the variable h over and letting e —)■ 0 yields 

the result. • 


This implies the following lemma: 
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Lemma 2.6 If(j){x,t) is either Ut, or (i G where u{x,t) is the 

dt\(j)\ + A“|0| < K\(j)\. 


solution of (1.1), then 


Additionally, we have 


dt\'tp\ + < |/"(m)#| + k| 


for'ijj{x,t) a solution to (2.5). 


Proof. Multiplying (2.3) by sgn(0) (respectively (2.5) by sgn('^)) and using Propo¬ 
sition 2^ yields the result. • 


3 The gradient decay estimate 

Throughout this section, all inequalities will be up to a constant independent of 
the solutions. We also weaken our requirements on the decay of the initial data. 
Specihcally, we insist that 


0 < Uo(x) < 


C 


1 -I- lx 


d+a+l ■ 


With the notations of Section 2, we would like to prove an exponential decay rate 


for the gradient (1.4) and second derivatives (1.5) of the form 


\4>{x,t)\ < u{Xif)e and \'4){x,t)\ < u{x,f)e (3.1) 

To do this, we will need to examine both representations for the evolution of the 


derivative ((2.3) and (2.4)). We have to distinguish three different ranges of |x|: 
long range (|x| » e^^), medium range (|x| ~ e'^*), and short range (|x| « 


-,\t\ 


The proof of (1.4) is stand-alone and presented fully. The proof of (1.5) is nearly 


identical (and we present it in parallel), but assumes that (1.4) holds a priori. 


3.1 Long Range: Kernel Estimates 

Lemma 3.1 For any fixed c > 0 and for any x in the range |x| > a 


solution 0 of (2.3) satisfies 


—CiXt 


(x, f)| < u{x,t)e 


Later, we will set the value of c based on the parameters of the problem. For now, c 


only affects the hnal constant that appears in estimate (3.1). As such, we suppress 
its notation in the proof of the lemma. We also draw special attention to the fact 
that, in the long range, we have an explicit exponent for the decay rate of Xa. This 
particular value will later facilitate a compatibility condition between the long and 
medium ranges. 
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Proof. We first prove the estimate for 0. Set g{u) = f{u) — ku. By Duhamel’s 
formula we have 


0(x,f)= dy^pa{x-y,t)uoiy)dy 


+ 



-y,t- s)dy^g{u{y, s))dyds. 


(3.2) 


Note that ||9xPa(-5is not integrable near f = 0, for a < 1/2. Hence the second 
term in the expression for 0 must be split into two pieces; one in a neighborhood of 
s = t (where jg essentially constant) and the rest where we may integrate by 

parts. In total, we obtain the following upper bound on 0(x, f): |0(a;,f)| < lo+Ii+hi 
where 

Jo = e'"* / uq{x- y) \dyPo,{y,t)\dy, 


-t-i 


h = 


/ \f{u{x-y,s))-Ku{x-y,s)\\dyPa{y,t-s)\dyds, 


'0 


l 2 = / \f{u{x-y,s))(j){x-y,s)-K(j){x-y,s)\paiy,t-s)dyds. 

Jt-l JR'i 

So we proceed by estimating each of the three integrals. We begin with Jq and apply 

r the 

id+a+l*^^’ 


the naive estimate for dxpa (2.8) and a concrete decay rate for the initial data uq. 

1 1 


Jo < e"‘P(f) 


1 + |a: - 1 + |9|‘ 


(3.3) 


where P{t) is a polynomial in t. The specihc form and degree of P are not needed 
since we will be getting a slight exponential decay in time, and that dominates any 


residual polynomial growth. To handle the main integral in (3.3), we split into two 


regions: {|?/| < |a;| /2} and {||/| > |a;| /2}. Observe that \x — y\ > |a;| — \y\ which, in 
the hrst case, means \x — y\ > |a;| /2. Thus 


-dy < 


dy 


<|.|/2 1 + |x - yr^^^ 1 + 1 + \y\ 


|d+Q( + l ^ l^jrf+Q + l ^ \ ^\d-\-OL-\-\ 

for X in the long range. For the second case, 

I' 1 1 , 

dy < 


< 


d+a+1 — I id+a+l 

tZ' 


-|>|x |/2 l + \x-y\ 


d+a+l 


1+ bl 


d+Q+l 


1 

f 

dy 


1 

X 

d+a+l 1 

1 + 

\d+a+l — 

\x-y\ 

X 

d+a+l 


Finally, observe that condition (2.1) and the size of |x| imply that u{x,t) 
^Kt 1 Since |x| > then have 

Jo < ^ u Ixp^ < 

which gives us our desired comparative decay for Jq. 














Continuing, we need to estimate the factors involving / and /'. Since / is smooth 
enough to have a partial Taylor expansion, we have that f{u) = /(0)+/'(0)m+O(m^). 


Therefore f{u) — nu = 0{v?), and this (along with (2.1)) gives one of the upper 
bounds we will use for Ji. 


pt-i 

h< 

Jo 


Q{t,s) 


1 + e-2«* 1 + |x - 2/1'^+“+^ 


dyds. 


(3.4) 


Here Q(t, s) is algebraic in t and s. As discussed earlier, Q{t, s) will not be integrable 
up to s = t, but it is at most polynomial in t on the interval 0 < s < t — 1. As with 
Jo, we handle the main integral by splitting into the same two regions: {||/| < |x| /2} 
and {|?/| > |x| /2}. In the first case, we get 


i<|x|/2l + + \x — y\ 


d+a+l 


dy 


< 


1 

f dy 

^Xds 

l^jd+a+l 


’1 |d+a+l 

\x\ 


(3.5) 


The last estimate follows from a change of variables. In the second case, we get 

r 1 1 


< 


i>|a:|/2 1 + 1 + k - 2/1 

f dy 


^2ks 


^ 2 k.s 


dy 


(3.6) 


\x 


2{d-\-a) 


1 + \x — y\ 


d+a+l 


\X 


2{d+a) ' 


Plugging the estimates from (3.5) and (3.6|) back into (3.4) shows that 


rt-l 


h < 


Q{t,s)e'^d-s) 


KS 

g d+Q. 

I d+Q+l 


^ 2 k,s 


+ 


1 2{d+0') 


ds < P{t) 


-.Kt 


^2K,t 


I d+a+l 


+ 


'0 \ |a;| Ixl ' ' / \ |x| ■ ' \x 

Remembering that u ~ e''* and |a:| > we see that 

^K.t 


2(d+a) 


h < uP{t) I 7^ + 


FI X 


d+a 


< uP{t) + 6-“'^*) , 


establishing the comparative exponential decay rate for R. Observe that our esti¬ 
mates here would not work if we only had |x| > e^*. 

We handle I 2 in an analogous manner. Note that f'{u) = /'(O) -|- 0{u), so that 


f'{u)(j) — Kcj) = 0{u(f)). But we also have, from Lemma 2.6 




xp)\ < l+Mol < Ct- 


\x 


d+a ’ 


and so |0(x,/:)| < Ctu{x,t), recalling that Theorem 2.2 implies u{x,t) > 

C'e'^*|a:|“F+«)_ This again implies a 0{tu^) upper bound which reduces the esti¬ 
mate on I 2 to 
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1 


l 2 < Q{t,s) 
Jt-1 


-dyds. 


, 1 + 1 + |x - 2 / 1 '^+“ 

Here, Q{t, s) is again algebraic in t and s, bnt can now be made integrable on [t — 1, t] 
dne to the correct time-homogeneity of pa{x, t). We estimate the convolntion integral 
exactly as we did for Ji and hnd that 


-t 


KS-t^ 

e 


^2ks 


h < l_QX‘>) [ ^ 1 


dXt 


^2nt 


+ 


|^|2(d+a) 


Continning as we did for Ji, we hnally get 

h < uP{t) = uP{t) (e-^ + . 

Finally, we examine the exponents to conclude that, for |a;| > we have 

\(f){x,t)\ < u{x,t)e~^'^^. Here, again, the method would not have worked if we 
merely required |a:| > 


Now we turn our attention to ^|J. Looking at (2.5), we employ the Duhamel formula 
to obtain 


i>{x, t) = J dypy^paix - y, t)u(i{y)dy 

+ [ [ -yp- s)dy.g{(f){y, s))dyds 

+ -y,t- s)f{u{y, s))((){y, sfdyds. 



(3.7) 


The comparative decay estimates follow analogously to the argument for 0. The 
hrst term has even better decay than /q since we placed two derivatives on p^- The 
second term is again split into two parts Ii and I 2 (for the time intervals [0,f — 1] 
and [t — l,t]). Since we have already concluded that 101 < in this range, the 

analysis for Ji proceeds identically. 

Now, f'{u)'ijj — ^0 = Qipitlj). The analysis of I 2 would also be identical to that of I 2 
if we knew that 101 = 0{u) (up to a time-dependent, not exponentially increasing 


factor). We prove this by appealing to Lemma 2.6 

|0(a;,t)| < *x \dxjd^^uo\ + - s) |/"(m) 00|(-, s)cis := Jq + Ji 

Jq is bounded by Cte'^^/\x\'^~^°' as before. Ji is estimated by the (now proven) 
comparative exponential bound for 0 (and 0) in the long range: 


Jl< 


< 

/d 

< 

F(f)e“^'^"* 

< 

F(f)e-2^"* {u 


l_,_g-2Ks|j^|2(d+a) l^^_y]d+< 


;dyds 


(3.8) 
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The third inequality was obtained through the same estimates as for I 2 . Notice 
how Ji inherits an exponential decay from the more rigorous analysis for 0; using a 


similar approach for the second term of (3.2) would have merely given us a compar¬ 
ative boundedness for 0 instead of decay. Hence Ji is dominated by Jq. Therefore 


ip = 0{tu), and we conclude that the second term of (3.7) satishes a comparative 
exponential decay rate in the long range. 


Lastly, observe that the third term of (3.7) is in fact dominated by Ji. As seen in 


(3.8), Ji also satishes a comparative exponential decay rate (with exponent —2 Aq!). 

This completes the mirrored argument, demonstrating that \ijj{x,t)\ < 

as long as |a;| > • 


3.2 Short Range: stability of the steady state u = 1 


The nature of / and condition (2.1) ensures that there exists a hxed (3 > 0 such 
that {|.^| < |/9} C {/'(m) < —(5o}- That is, we have an exponentially growing ball 
wherein u is sufficiently close to 1 that f'{u) is negative by at least a hxed amount. 
Let Ho(f) = {|.^| < |0}. In this region, we use a simple version of the maximum 
principle. Assume the positive maximum value of 0 occurs at time t at an interior 
point X E Qo{r) (negative minimum values are handled similarly). Then 

0t0(x, t) + A“0(:r, i) = f'{u{x, f))0(x, f). 

However, we have A“0(a;, f) > 0. Therefore, 


dt(l){x,i) < -5o(j){x,f}. 

We conclude that either 0 is decaying exponentially (with exponent at least —5o) on 
or the size of 0 on Qo{t) is controlled by the size of 0 outside Ho(f). Equivalently, 
V is eventually controlled by a hxed constant times whenever |^| < |0 so 

long as V obeys the same bound for (say) |^| > 0, which will be established in the 
next two sections. Certainly, the previous argument shows that this is the case in 
the long range. 


To prove the analogous statement for 0, we will assume that (1.4) holds with some 
uniform 5 > 0 in all ranges for 0 (the proof of which will, as stated, be stand-alone). 
Since u is bounded by 1, we have that 


dti>{x, t) + A"0(a:, t) < —So^jJ{x, t) -)- Ce 


-2St 


for all X G Assume the positive maximum value of 0 occurs at time t at an 

interior point x G Ho(^)- Then 

dt'ip{x,t) < —Soip{x,i) + 

again establishing that either 0 is decaying exponentially on M'^, or the size of 0 on 
f2o(^) is controlled by the size of 0 outside Ho(i^)- As before, negative minima are 
treated similarly. 
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3.3 Medium range: selective comparison principle 


For the intermediate range, we appeal to Theorem (2.4). Let t2i(t) = {/3 < |^| < 
ce^"*}, for c > 0 to be determined later. We need to hnd a w such that 0) > 

\v{^, 0)1, and w > |n| outside of for all f > 0, and such that 

dtw — ■ Vw + — Xw — f{u)w > 0 (3.9) 

pointwise on Vli{t). It is important that w be dehned on all of because A“ is 
nonlocal; while it can be dehned on compact sets, it becomes a different operator 
than the one in (2.4). 


We must therefore provide a candidate supersolution w. For ^ G and 

0<i> — d — a<l, dehne the family of functions Wy by 

Mi. t) = ^ 

where Cy,WQ > 0 will be determined later. Observe that each such Wi, is radial. To 
simplify notation, let y = |^|. For hxed u and t, , Wu is the solution of the following 
ODE in y (on Di(f)): 

y-^Mv, t) + i^Mv. t) = (3.11) 

dy 

Note that = Cu/iv — d + a). Cy will be determined below. Technically, Cy 
will also depend on wq, but the dependence will be linear. More importantly, this 
dependence is not circular: wq is the initial condition aX y = (3, and this may be 
hxed before solving the ODE. wq itself will depend on the size of 0). 


Our supersolution must be positive, yet Wy is eventually negative for sufficiently 
large y. We now choose c such that tCi, > 0 on Di(f); since z/ < d + a + 1, such a 
choice is always possible. Outside of Di(f), we extend Wy to be smooth, bounded, 
positive, radial, decreasing in y, and This can obviously be done such 

that the C^, and norms of w are uniformly bounded for all f > 0. They do, 
however, depend on u and wq. 


We now dehne our candidate supersolution as w{^,t) = e^^Wy{\^\ ,t) with rj and 
u hxed; we will determine their values shortly. Observe that ^ ■ Vw = ye'^^dyW. To 
prove (3.9), it therefore suffices to show that 




hWy 


yw + Xydywiyj, f) + (A + /'(m(^, i))) t) < -e Mw{y, t) \ 


on Di(t). Examining (3.10), we see that dtWy > 0 on our domain, so we get a 
stronger inequality if we ignore the hrst term above; proving the stronger inequality 
will be sufficient. Since f'{u) < k = X{d + a) and w > 0, we can eliminate the ^ 
dependence entirely and satisfy the above inequality provided that 
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Xydyw{y, t) + X{d + a + 1 - j)w{y, t) < -e \A°‘w{y, t) \. 

Treating the diffusion as a perturbation term, we need to bound the size of 
\A°‘w{y,t)\. Our choice of extension for Wu onto implies (see [3], and estimate 

(8) of m-- 


\A'^w{y,t)\=e^^\A<^w,{y,t)\<e^^ 


q 


on r2i(t). Our assumptions on the uniform-in-time bounds of the norms of w,, ensure 
that the constant q can be made time-independent (recall that u is hxed). It is 
also linear in wq, as mentioned before. Thus, proving (|3.9) reduces to showing that 




ydyw{y, t) + {d + a + l- -^)w{y, t) = 


_ _ ^rjt—Xat 


d-\-a ' 


y 


with q > 0 now determined. But this is precisely (3.11), multiplied by with 


i' = d + a + l — j. Hence, (3.9) holds on Oi(t). Clearly, we need 0 < rj < X. We 


then choose Wq sufficiently large that w (that is, the function extended to all of ^ 
dominates Ini at t = 0. 


In order to use Theorem (2.4), we need to show that tn — |n| > 0 outside of It 

is here that we needed our free parameters u and y, naively, we would like to take 
T] = 0 and u = d + a + 1, SiS this would give the best bounds in the medium range. 
However, the upper bounds for cj) established in the previous two sections are too 
weak to ensure that w — |n| > 0 outside of for all t with that naive choice. 

In the long range, we only have |n(.^,t)| < the short range, we 

have that |n(.^,t)| < eA~^o)t pj^ovided the same can be said of v globally. 


The added growth rate (in the form of alows us to ensure that w > |n| in 
the long range. Since both w and the upper bound for v are radially decreasing 
and 0(1^1”^'^"'"“^), it suffices to check the inequality for |^| = At the boundary 

between medium range and long range, w ~ ^ ^-(d+a)\at+M-\at^ 

So it suffices to have rj — Xa{u — d — a) = y — Aa(l — j) > A(1 — a), or y > 

This leaves enough room for y to be smaller than A. Fixing constants (and looking 
at the equation after a transient period to alow the long range decay estimates to 
come into effect) we have that w > |u| for |,^| > ce''‘“*. 


For the second derivatives, we will similarly employ a supersolution of the form 
= e^*q(|^|,t) with nearly the same q as in (3.10); specifically, we require 
a larger constant q in equation (3.11). Essentially, we need W to satisfy 


- dtW + A^ ■ VIE + 2A1F + f{u)W < -e-^“*|A"hF| - (3.12) 

once again appealing to the full estimate for <p. As before, |A“hF(^,t)| < 
qe‘“*|^|“^‘^’''"^. But the last term is dominated by (7e2A-25t|^|-2(d+o)_ jp gg_ 
sume that ( > 2X — 26, the diffusive term is then a stronger perturbation than the 
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inhomogeneity, and so the latter can be absorbed into the bounds for the former. 
Our candidate supersolution W will solve (3.12) provided that 


C C',y 

ydyW^{y, t) + {d + a + 2- pw^{y, ^ " 


d+a ’ 


This is just (3.11) with i/ = d + a + 2 + (/X, so the equation is satisfied. We want to 
take C as small as possible while still observing the requirement that 0 < u—d—a < 1 
(this is necessary to ensure W remains positive on fli{t) with uniform bounds). Thus 
our arguments impose the constraint that 2X — 26 < ( < 2X. 


For the final constraint, W must dominate V (that is, in exponential co¬ 

ordinates) at the interface between medium and long ranges. As before, we have 
that W ~ g-J^Aat+C* 1^1 < ^-{d+a)\at+ 2 Xt-\at |^| = We therefore 

require that ( — Aa(z/ — d — a) = ( — Xa{2 — ^) > A(2 — a), or (^ > Aj^. For a > 0, 
this will always overlap with the admissible range of ( determined above; and, as 
expected, the minimum value for ( decreases as we get better diffusion. 


This establishes, as before, that IF > |1A| in the long range and that W satisfies the 
comparison inequality on Qi{t). To satisfy all the conditions for the selective com¬ 
parison principle, we must show that the supersolutions for v and V also dominate 
in the short range, f2o(^)- 


3.4 The estimate for (J) and tjj 


Proof of Theorem |1.2[ Recall that = {|^| < §/?} and Qi{t) = {f3 < |^| < 
Section 3.2 established that, in f2o(t), v < or else v is larger 

outside of this regime; which, by Lemma 3.1, in fact means v is larger in r2i(f); 
outside of h2o(^) U r2i(f) (a compact set), the comparative exponential bound holds 
unconditionally. 

Fix r] = max (A — 5o/2) and Wq so large that the conditional bound in the 
short range is stronger than the requirement that w > v in the same regime; that 
is, for all t > 0 (( |3.10[ ) shows that such a choice is 

always possible). We can also ensure that w > v at t = 0. 


We then argue by contradiction. Let f > 0 be the critical time such that 
for all t < f and ^ G f2o(^) \ and such that the same 

upper bound fails for a sequence t t"*"- This therefore implies a sequence 
{^n} e ^o{tn) \ such that v{^n,tn) > w(^n,^n) > C(no)e^^"'^°^*. By the 

conditional bound of Section 3.2, there is another sequence G f2i(tn) \ f^o(^n) 
such that > inf^gfjQ(i^)\n^(t^)t(;(,^, f„). By compactness and the continuity of 

V, we extract a limit point ^ with |,^| > |/1 and n(^, ^ > inf|^|<^w(,^, f). 


However, for t < t, the hypothesis of Theorem (2.4) was by assumption valid, 
and therefore v < w up to t. Moreover, (3.10) immediately shows that there is a 
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continuous, always positive e{t) such that sup|g|>3^w(.^, t) < inf|^|<^w((^, t) — e{t). 
The continuity of v allows us to extend the inequality to t. Thus 


u(^, t) > inf|g|<^'u;(^, i) > sup|^|>|^'«;(^, i) + e{i) > w{^, i) + e{i) > v{^, i)+e{t), 

which is evidently a contradiction. Therefore, w > v (and analogously w > —v) 
in \ The hypothesis of Theorem (2.4) is always satished, and so 

> |n(.^,f)| for all ^ and t. Precisely the same argument also applies to V 
with ( = max (Aj^, 2A — 2S, 2A — (5o) • 


Putting everything together, we have the following three upper bounds (written 
in terms of 0). Recall that u ~ in the medium range. 

\(j){x,t)\ < C 2 u{x,t)e~^°'^ for |x| > 

\(j){x,t)\ < Ciu{x,t)e^^~^^ < Ciu{x,t)e~^^^ for cqc'^* < |a;| < 

\(j){x,t)\ < Coe~^°^ < Cou{x,t)e~^°^ for |x| < cqc^^ 

We also get the corresponding three bounds for ■0 (recall that ip = 

\ij{x,t)\ < C'^M(x,t)e-^“' for |a:| > 

\'ip{x,t)\ < C[u{x,t)e‘^^~^^^ |^|-2+C/^ ^ C[u{x,t)e~^^^ for cqc'^* < |a;| < 

\'ilj{x,t)\ < Coe~^°^ < Cou{x,t)e~^°^ for |x| < cqc^* 


Thus proving (3.1) and Theorem (1.2). 


Remark: Theorem (1.2) can be extended to a comparative exponential decay for all 


derivatives of u. The dual proof presented above is inductive, essentially repeating 
the same arguments for ip while assuming the full result for 0. This can be reworked 
into a full induction showing the much stronger result that 


\'V^u{x,t)\ < Cu{x,t)e 


Sift 


for all k. The next section will, however, only require (3.1). 


4 Bound on the fractional Laplacian 

We hnally come to the estimate telling us that, at large times, the fractional Lapla¬ 
cian of the solution u is small compared to the size of u. 

Lemma 4.1 There exist C > 0 and > 0 such that 

\A°‘u{x,t)\ < Cu{x,t)e~^'^. (4.1) 


15 







Proof. We will examine the explicit integral representation for A°'u{x,t) and ini¬ 
tially with |x| > ce^*. Up to a constant, we have 


A°‘u{x,t) = 


dy + 


\x—y\<ee'^^^ 

h + h 


' \x—y\>£eA^^ 


dy 


u{x,t) - u{y,t) 


\x 


y\ 


d+a 


(4.2) 


where 7 G (0,1]. For the inner piece Ji, we make use of (3.1) and a Taylor expansion: 
u{x,t) — u{y,t) = —Vu{x,t) ■ {y — x) + O {\V‘^u{cy + (1 — c)x,t)\\x — j/p) for some 
c G (0,1). Keeping in mind that 


x-y 


Ja \x-y\^+^ 

on any annulus A centered at x, we then have 
I V^u(ca; + (1 — c)y, t)| 


dy = 0 


\h\ < 


-dy < e 


-5t 


u{cx + (1 — c)y, t) 


' |x —\x y\ 


d+a —2 


dy. 


I I I d+Q— 2 

J \x—y\<£e'y^'^ \X y\ 

Now, in this regime, u{x, t) ~ e''* and |cx + (1 — c)y\ > (1 — e) |x| because 

we assumed \x — y\ < < e\x\^ < e\x\. Therefore, 


lUl < (1 






-^d—a —2 


dr. 


We then have (up to a constant depending on e) 

|/i| < < u(x,t)e-'’‘l'‘. 


(4.3) 


To obtain the last inequality we must choose 7 appropriately. If (2 — a)A < 5/2, we 
may take 7 = 1 . Otherwise, we take 7 = 5/(2A(2 — a)) < 1. Either way, we obtain 
our final estimate for Ji. 


For the outer piece J 2 , we merely use the positivity of u\ {u{x, t) — u{y, t)) < u{x, t). 


Ihl < u{x,t) 


dy 


\x—y\>ee''!^^ \x y\ 


d+a 


= u{x, t) 


coo ^d 1 ^^ 


e7At r 


d+a 




The above inequality and (4.3) show that |A"u(x,f)| < for some fixed 

positive 5' and |x| > ce'*'*. The exact value (in terms of 5, a, and A) requires an 
optimization argument for 7 . 


The case where |x| < ce^^, i.e. the short range, is much simpler. In this range, we 
have that 0 < c < u{x,t) < 1- And, moreover, we know from Theorem 1.2 that 
llv+IU 00 < Ce Using our Taylor expansion again, we can interpolate between 
these two to obtain: 


u{x,t) - u{y,t) 


\x - y\ 


d+a 


dy 


< 


|V+(.,t)i|Locd?/ 


2dy 


'\x-y\<R |x — y\ 


d+a—2 


'\x-y\'>R \x — y\ 


d+a — 




The last inequality required us to optimize R to || V^m(., and use our uniform 

estimate. Hence \A°‘u{x,t)\ < Ce~^^ < Cc~^u{x,t)e~^ ^ for all x. • 
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5 Asymptotic symmetrization 


We have used ( |3.1[ ) to establish a similar comparative decay on the asymptotic size 
of A“m. We will now look at the explicit (integral) representation for the solution 
u{x, t) and obtain a better approximate asymptotic behavior than what may be 


inferred from (2.1). We will then use this to prove Theorem 1.3 


5.1 Behavior of u at infinity for t = 1 

The second ingredient in the symmetrization proof is the fact that u{x, 1) decays 
exactly as for large x. This is detailed in the following lemma. 

Lemma 5.1 For all t > 0, there is kt> and 5 > 0 such that 


u{x, t) = 


I ('i+«) 


\X 


+ 0 (- 


^ I d~\~cx~\~S 


). 


(6.1) 


Proof. Once again it is a matter of applying Duhamel’s formula to u: 
^Kt /* ^r 


u{x, t) = 





:P- 


x-y 


{t - sYi°^^°‘^{t - syF 


Id t / Jq 


W{y,s)dy 




We have, by Proposition 2.3 


because uq is exponentially decreasing. As for the second term, it is sufficient to 
estimate D{s,t,x,y), the integration in time being only between 0 and t. Let us 
decompose the integral 


D{s,t,x,y)= +/ ;=/i + / 2 . 

J\x-y\<\x\/2 J\x-y\>\x\/2 


We omit, in order to keep the notations light, the dependence in s, t, x, y. As for 
Ji, let us notice the bound u{x, t) < -j—simply because we have dtU + A^^u < nu. 
So, we have 


h < 
< 


Ct 


\x\Ad+o^){t-sYI- y|x-,|<N/2 

Ct II II 


Pc 


I I 2 (d+Q:) 


X-y 

{t-syc 


)dy 
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As for / 2 , let us use estimate (2.9) in Proposition 2.3 and write I 2 < I 21 + ^*-^ 22 , 
where 

u^{y,s) , . Ct 


I 22 = {t- s) 


and 


J|a:-i/|>|x|/2 

I 21 = {t ~ s) [ 


\x — y|'^+"+'5 


dy < 


I ^ I d-\-0!.-\-S ^ 


'|x-y|>|x |/2 


u^{.y,s) 
\x — y\d+a 


dy. 


Let us split I 21 as 


L 21 — / + / — L 211 + / 212 , 

J\x-y\>\x\/2,\y\<e\x\ J\x-y\>\x\/2,\y\>e\x\ 


c, 


t.e 


where e > 0 is small. We have /212 < , ,r,r ,, n 
W e use the fact that, for \y\ < e\x\ we have 


\x-y\ = |x|(l + 0 

in other words we have, setting 

:= {|x -y\> \x\/2, \y\ < e\x\} 
and using once again estimate (2.9|): 


, so it remains to study the last term. 


\y\ 


I 211 — 


Ca{t - S) 


./c, 


u‘^{y,s)dy + 0 


+0 


{t - / 2 


t — s 


\^^\^d+a+l 


\y\u^{y, s)dy 


I ^ I d~\~cx.~\~S 


u {y,s)dy 


'c^ 


Set l{t) = / u^{y,t) and note that / \y\u^{y,t)dy < + cxd . Hence there is > 0 


such that 


_ ^a{t s)ls 1 

^211 - |rf+„+ 5 h 


X 


Gathering everything yields (5.1). 


5.2 Analysis of long-time level sets 


The proof of Theorem 
Proof of Theorem 


1.3| now reduces to the explicit resolution of a simple ODE. 


1.3, For all A G [0,1), there exists mx > 0 such that 


Vn G [0, A], mxu < f{u)/K < u. 


The combination of this remark and Lemma 4.1 reduces the full problem (1.1) to 


dtu = (1 + 0(e '’*))/(«), u{x, 1) = + 0(—as |x| ^ +cx), (5.2) 

Jb 


\X 
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keeping in mind that 0(e might depend on A. Set 


G(u) 


I 


du 

wr 


as M —)■ 0 we have 


G{u) = +go + 0{u), 


K 


where 5^0 is a positive constant. Finding the A-level set of u for large t amounts to 
integrating (5.2) with x large, which yields 

r dtu 




'1 f{u) 


in other words 


Logn + 0{u) — Log 


ki 


\x 


d+a 


+ 0 { 


^ I d-\-(y.-\-S' 


) ] — Kt qoo{uo) + 0{e 


Specializing u{x,t) = A yields 


exp 0( 


\x 


d+a ' 


ki 


\x 


\d-\-a. 


+ 0 (- 






for some possibly different constant qx^ooiuo). This, after elementary computations. 


yields (1.6). 


Let us point out again that we have proved, in quite a strong sense, that the large 


time dynamics of the KPP problem (1.1) is the same as that of the ODE ii = /(m). 


the fractional diffusion being there only, but this is important to set the initial da¬ 
tum right. This fact was already noticed in [6] and, in an even more evident fashion. 


m 
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